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Abstract. Global star formation is the key to understanding galaxy disk formation. This in 
turn depends on gravitational instability of disks and continuing gas accretion as well as minor 
merging. A key component is feedback from supernovae. Primary observational constraints on 
disk galaxy formation and evolution include the Schmidt-Kennicutt law, the Tully-Fisher relation 
and the galaxy luminosity function. I will review how theory confronts phenomenology, and 
discuss future prospects for refining our understanding of disk formation. 
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1. Introduction 



Or 

: 

, There are three interrelated empirical scaling functions that relate all disk galaxies and 
are fundamental to understanding their formation. These are the Tully-Fisher law, the 
Schmidt-Kennicutt law, and the galaxy luminosity function. If we can ever understand 
these, we will have come a long way towards understanding disk galaxy formation. I will 



^ . 

review some of the issues that arise when theory confronts phenomenology. Star forma- 
tion and mass assembly are the primary issues to be understood. Although intimately 
connected with star formation, I will mostly avoid chemical evolution constraints in this 
discussion, in part because of their strong sensitivity to adoption of an IMF and because 
l — they are covered elsewhere in this conference. 

^ . 2. Tully-Fisher relation 

t-H , The Tully-Fisher relation is a long known correlation between disk luminosity and max- 
imum rotational velocity. The principal current issues in understanding the Tully-Fisher 
(TF) relation involve angular momentum loss and adiabatic contraction. The I-band 
0*N Tully-Fisher relation, preferred to the B-band data because of the lower scatter, provides 
a challenge for understanding current epoch disks. The contribution of turbulence and 
rotation to line widths reduces TF scatter and is measured to z ~ 1. No significant evo- 
lution in the TF relation is seen (Kassin et al. 2007). Even galaxies as far back as z ~ 2 
(Bouche et al. 2007) remain on the TF law as measured locally. 

Loss of angular momentum of the baryons to the dark halo was discovered in numer- 
ical simulations (Steinmetz & Navarro 1999). Strong SN feedback can reduce angular 
momentum loss (Governato et al. 2007). The observed normalisation can be obtained in 
the inner disk. This success comes at a price, however. One concern is that of the impli- 
cations for the galaxy luminosity function at low luminosities. Another is disk thickness. 
The simulated disks are thicker than observed disks although the numerical resolution is 
inadequate for a real confrontation with data. 

The radial structure is a more pressing issue. Lack of evolution is observed in disk sizes 
at z < 1, in contrast to the remarkable compactness found at z ~ 2 — 3 (Buitrago et 
al. 2008). Adiabatic contraction of the dark halo due to the self-gravity of the collapsing 
baryons helps account for this by modifying the effective size evolution at early epochs 
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(Somerville et al. 2008). However the response of the dark halo to baryon dissipation by 
adiabatic contraction means that once the rotation curves are properly corrected, the 
asymptotic rotation velocity fails to match the observed TF zero point (Dutton et al 08). 
An equivalent symptom of the TF issue is that the predicted mass-to-light ratio is too 
high. 

3. The Schmidt-Kennicutt law 

The disk star formation rate is coupled to the gas surface density by the empirical 
relation 

SFR = 0.02(GAS SURFACE DENSITY) / td yn ■ 
Remarkably, this relation fits quiescent as well as star-bursting galaxies. The individual 
star- forming complexes in M51 also lie along this line (Kennicutt et al. 2008). This is more 
of a mystery and suggests that local processes control the efficiency of star formation. 

We may regard the Schidt-Kennicutt (SK) law as specifying the disk mode of star 
formation. It is motivated globally by the gravitational instability of cold disks. Feedback 
from supernovae is crucial to avoid too high a star formation rate. Momentum balance 
yields 

SFE = a gas v cool m* SN /Egl ial « 0.02. 

The global injection of momentum is balanced by local dissipation. Here m* SN is the mass 
in forming stars required to produce a supernova of Type II, approximately 150 M for 
a Chabrier IMF, and E^ tial « 10 51 ergs is the initial kinetic energy in the supernova 
explosion. 

In this way, a near steady state of cloud motions is obtained at the disk interstellar 
cloud velocity dispersion of a w lOkm/s. This corresponds to the peak in the diffuse 
ISM cooling curve. The local control of momentum input and dissipation is to be found 
in exploring cloud disruption once massive stars form, develop HII regions and explode. 
The explanation is derived from the local physics of cloud efficiency which controls the 
efficiency of star formation, that is, the fraction of gas converted into stars within a 
specified time-scale. The relevant local time-scale is molecular cloud lifetime. The com- 
bination [(SFR x dynamical time) /gas density is approximately constant as a function 
of gas density, and amounts to about 0.02 (Krumholz & Tan 2007). 

There are several unexpected developments for the Schmidt-Kennicutt (SK) law. Star- 
forming complexes in M51 individually lie on the SK law, albeit with increased dispersion 
and a slight steepening of the slope (Kennicutt et al. 2007). Extreme starbursts fit the SK 
law both in the nearby universe (Sajina et al. 2006) and to z ~ 2 for SMGs (Bouche et al 
2007). The SK law is dominated by the molecular content of the galaxies. This is not too 
surprising as stars form in molecular clouds. The total molecular mass measured via CO 
amounts to ~10 % of the HI in LMC, M33. In the MWG, the atomic and molecular gas 
masses are comparable, and many ultraluminous infrared-emitting galaxies are molecule- 
dominated (Gao & Solomon 2004). Use of molecular indicators that trace the densest 
gas such as HCN reveal a linear law whereas use of CO, which traces the least dense 
molecular gas, prefers the SK law (Komugi et al. 2005). The change of slope is most 
simply understood in terms of the excitation thresholds (Krumholz & Thompson 2007), 
with the dominant mass in molecular gas, traced by CO, representing the physics of the 
self-gravitating cold disk instabilities and cloud collisions that drive the SK law. 

A more quantitative model of ISM turbulence comprises cloud collisions, with stirring 
by gravitational instabilities and by SNe momentum injection. Let a typical cloud have 
pressure p c i and surface density S c ;. We expect star- forming clouds to be marginally 
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self-gravitating, so that p c i — kGY? c1 . They are also marginally confined by ambient gas 
pressure, p g = p g <J 2 . Clouds certainly form this way, and this state may be maintained 
if the cloud covering factor is of order unity. This guarantees that cloud collisions occur 
on a local dynamical time-scale. We have 

p g = nGY, g Y, tot = pd = 7rGE c / 2 , (3.1) 

from which S c j = (T, tot 'E g ) 1 ' 2 . Then the covering factor f2 c ; of clouds in the disk is 

1/2 

directly inferred to be f2 c / = S 9 /E c / = (X g /£ tot ) ' . The cloud collision time-scale is 
t co u — (E c /i?)/(Sgcr g ), where the scale height H^ 1 = (nGY.tot) / cr 2 , and a g is the cloud 
velocity dispersion. The collision time Cclll cllsO be expressed clS t co ll — Q>cl ^cross 

with 

1 /2 

t cross = H/o-g, which becomes t co u = (£ tot /£ s ) (H/a g ). 

Assume the disk star formation rate is self-regulated by supernova feedback which 
drives the cloud velocity dispersion. Momentum balance gives 

£*(E S N/(m* SN v c )) = f c Zg(jg/t co u. (3.2) 

Here f c is the cloud volume filling factor, which can be expressed in terms of porosity Q 
as f c = e~®. Also v c is the velocity at the onset of strong cooling. Canonical numbers 
used throughout are m* SN — 150M© and v c = 400 kms -1 . 
We can rewrite this expression as 

P* = tSNfc\/Gp g pg (3.3) 

with 6sn — {m* SN v c a g )E S N~ 1 (Pg/Pci) 1 ^ 2 ■ For disks, we obtain 

S* = M(irG) (S tot ) V2 ((msNVc)/E S N)}[p g /pci}^ g 3/2 

= e SN f c ^T g (R/H)^ gas n, 

where Y>tot — S s + S*. Here the disk gas fraction is f g ~ 0.1, and we use the disk 
radius-to-scale- height relation H/R — (o- g /v r ) 2 for a disk rotating at v r with Q 2 = 
GY>tot I R- Remarkably, although the preceding formula ignores the multi-phase nature of 
the interstellar medium and the possibility of gas outflows (see below), one nevertheless 
manages to fit the Schmidt-Kennicutt relation (Silk & Norman 2008). Mergers, which 
certainly drive starbursts, will increase the cloud velocity dispersion. However there is self- 
regulation between star formation efficiency e$N fx cr and and disk scale-height H oc a 2 . 
This accounts for the robustness of the SK law fit. 

There is one important addition to the above discussion. The disk heats up dynamically 
unless a supply of cold gas is provided. Cold gas accretion is needed in order to maintain 
global disk instability and continuing star formation. Major mergers cannot provide the 
gas source, as they would be too destructive for thin disks. Minor mergers provide a 
motivated source. Observations provide strong justification for an extensive reservoir of 
neutral gas, as seen in the M33 group (Grossi ct al 2008). Many nearby spirals are found 
to be embedded within extensive envelopes of HI. In some cases, such as that of NGC 
6946, there is direct evidence of angular momentum loss as the HI accretes into the star- 
forming inner galaxy as well as of extensive fountain activity (Boomsma et al. 2008). In 
all cases, the low star formation efficiency is plausibly due to SNe feedback, although this 
may not be the only mechanism at work. For example, magnetic pressure may also play 
a role in globally controlling gas flows via the Parker instability. 

3.1. The rate of star formation 

In reality, we find a multiphase interstellar medium which consists of hot phase (~ 10 6 
K gas) as well as atomic and molecular gas. The heating and cooling of this gas, as 
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well as mass transfer between the different phases, controls the rate of star formation. 
A simple porosity description of supernova feedback in a multiphase ISM (Silk 2001) 
provides an expression for the star formation rate in which porosity-driven turbulence 
is the controlling factor: p* = Qm5 A r(47r/3i?^t a ) _1 . The shell evolution is generally 
described by (Cioffi, Bcrtschinger and McKee 1988) t = toE^rig 4/7 (v c /v) w / 7 and 
R = RqEI' 7 n^' 7 {t/t Q f/ w , where v a = 413km/s, R = 14pc and t = 1.3 x 10 4 yr. 
Here v c = Al'iEri n„ 

/4 A 3/8 k 

m s where the cooling time-scale within a SN-driven shell 
moving at velocity v c is t c — v c /\p, A -1 = 3TOp^ 2 fc 1 / 2 T 1 / 2 /A// and A e //(T) is the 
effective cooling rate ( £ i -1 / 2 over the relevant temperature range). 

The SNR expansion is limited by the ambient turbulent pressure to be p g o~ g , and we 
identify the ambient turbulent velocity dispersion v a with a g . We obtain 

P, = Qc^m* SN nf/ 7 E-^ 14 (a g /vo)^ 7 
where Co = ^ L i?gto- Rearranging, we have 

P* = Q^Gp~ gPg (a g /a ftd ) 19 / 7 , (3.4) 

where 

« 20n 3 - 1 /i4 m ^37 oi? O i 47 km/Sj 

where m* SN is normalized to 100 M Q and E$n to 10 51 ergs. We rewrite the star formation 
rate as M* = 6QM g /td, where eg = Q(a g / 'a fid) 2 ' 1 ' fg^ 2 '■ Now feedback physics compels to 
identify eg with €$n- It follows that Q oc cr -1 - 7 . Hence porosity decreases with increasing 
turbulence. 

Simulations of star formation with SN feedback in a multiphase ISM confirm that 
porosity provides a good description of feedback. The analytical formula fits the numerical 
simulations performed at high enough resolution to follow the motions of OB stars prior 
to explosion (Slyz et al. 2006). The star formation rate prescription has been tested in 
numerical simulations of a merger. In the case of the Mice (NGC 4676 a,b), a pure density 
law fails to account for the extended nature of star formation. A turbulence prescription 
for the SK law gives a better fit to the observed distribution of young stars (Barnes 2004) . 

3.2. Application to starbursts 

Starbursts are characterised by high turbulent velocities and strong concentrations of 
molecular gas. Locally, the gas density where stars form is characterised by molecular 
clouds. It is presumably the cloud concentration that is enhanced by the dynamics of 
tidal interactions, which includes both stirring by bars and mergers. Now the porosity 
ansatz, if porosity self-regulates, leads to an expression for the rate of star formation: 

SFR = POROSITY x (TURBULENT PRESSURE) 136 . 

We have already noted that if turbulent velocities are high, the porosity must be low. 
This conspiracy, together with that between gas disk scale-height and star formation 
efficiency parameter, helps explain why even extreme merger-induced starbursts stay on 
the Schmidt-Kennicutt law. 

It turns out that the molecular gas fraction is regulated by the turbulent pressure (Blitz 
and Rosolowsky 2006). The molecular fraction is found to be approximately proportional 
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to the pressure, so that the star formation rate empirically is 

/ \ 0.92 

S 5F « = 0.1e£ 9 (^J Mopc^Gyr- 1 . 

Here e is an empirical SK law fit parameter. The molecular fraction is (p m oi/po) ' 92±0 ' W , 
where po/k — 4.3 ± 0.6 x 10 4 cm 3 K. Since S s p^ ol , the inferred star formation rate 
dependence on pressure is similar to that obtained at constant porosity. Disk simula- 
tions that include a prescription for star formation in molecular clouds (Robertson and 
Kravtsov 2008) can reproduce the slope and dispersion of the SK law.. 

It is plausible that porosity should self-regulate. Low Q is associated with enhanced 
turbulent pressure. The supernova remnant expansion is blocked, further enhancing the 
pressure and squeezing the clouds. This triggers more star formation and supernovae 
that result in blow-out of fountains from the disk. Blow-out drives high Q. If Q is high, 
we obtain lower pressure and infall is allowed. The infall in turn reduces Q, and drives 
further accretion that eventually results in star formation and supernova remnants. This 
in turn enhances Q. 

Another reason that Q self-regulates is that it is insensitive to the strength of the 
turbulence. The volume filling fraction of cold gas depends only logarithmically on the 
Mach number (Wada and Norman 2007). 

The resulting star formation is not expected to be monotonic in a porosity-driven 
prescription. The time-scale for variation is of order the life-time of a molecular cloud 
complex to disruption by SN feedback, of order 10 7 — 10 8 yr. A local burst-like behaviour 
is indeed seen in the solar vicinity (Rocha-Pinto et al. 2000). The inner star- forming 
regions of nearby galaxies also reveal non-monotonic behaviour. For example, M100 has 
undergone a sequence of bursts (Allard et al. 2006). 

Supernova-driven outflows lead to entrainment and loading of the cold ISM into the 
fountain and/or wind. The degree of loading depends on the ability to resolve Kelvin- 
Hclmholtz instabilities. This is completely inadequate even in state-of-the-art simula- 
tions. Empirically, the load factor is a few. as measured in NGC 1569 (Martin et al. 
2002). 

The fate of the outflowing gas is uncertain. The outflow rate is proportional to Q times 
the star formation rate times the load factor. Since Q is suppressed as er -1 ' 7 , it follows 
that outflows are suppressed in massive galaxies. One still expect fountains to recirculate 
gas from disk to halo and back to disk, as the halo gas cools. 

Star formation is observed to continue below the Toomre threshold for disk gravita- 
tional instabilities. This is naturally explained in the cloud collision model where even the 
outer regions are dynamically populated with orbiting disk clouds. The star formation 
rate is most likely modulated by the molecular gas fraction which itself is controlled by 
the local UV radiation field (Schaye 2004). 



3.3. Some current issues in understanding the Schmidt- Kennicutt law 

We still need to attain a deeper understanding of the universality of the SK law. It 
prevails both locally as well as globally, and at high as well as at low redshift. It applies 
in regions of extremely high star formation, such as extreme starbursts, and in regions 
of very low average star formation rate, albeit somewhat suppressed, such as the outer 
regions of disks and DLAs. 
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4. Galaxy luminosity function 

The final diagnostic, and relic, of disk formation that I will discuss is the galaxy lu- 
minosity function. The case for outflows has been made in the context of the cold dark 
model, as a means of suppressing the numbers of dwarfs (Dekel and Silk 1986) and for 
driving early chemical enrichment of the intergalactic medium as inferred from the LBG 
"missing metals" problem (Pcttini et al 1999) and wind-driven enrichment simulations 
(Cen & Ostriker 1999, Oppenheimer & Dave 2006). Early numerical simulations of indi- 
vidual winds (Mac Low & Ferrara 1999) argued for incomplete mixing of supernova ejecta, 
and this problem has not yet been resolved via multiphase simulations. Suppression of 
winds in massive galaxy potential wells was demonstrated by Springel and Hernquist 
(2003). Ejection of baryons may be accomplished either in the assembly stage via dwarf 
winds, as is relevant here, or from massive spheroids via AGN-triggered outflows. The 
latter is a separate topic that itself merits an entire review. 

An alternative, and complementary, means of dwarf suppression appeals to rcionisation 
inducing mass loss from low mass dwarfs, with escape velocity £ 20 km/s. This has been 
implemented in simulations but fails to account for the shape of the observed Milky Way 
luminosity function of dwarf satellites within 280 kpc (Koposov et al. 2007). The observed 
effective power-law index is a = 0.25 from My = —2.5 to —18. Either there is a poor fit 
to the observed luminosity function (Somerville 2002), or the predicted survivors have 
far too high a central surface brightness (Benson et al. 2002). Tidal stripping cannot 
account for the discrepancy (Penarrubia, Navarro & McConnachie 2008). 

There is an equally serious problem for giant galaxies. The problem here is that ac- 
cretion of gas continues over a Hubble time, and the galaxies continue to grow. The net 
result is too many and too blue massive galaxies. Supernova feedback cannot resolve 
these problems, but AGN/quasar feedback is considered to provide the needed panacea 
by heating the infalling gas and quenching star formation. The relevant feedback is that 
associated with supermassive black hole growth, and is destined to yield an explanation 
of the Magorrian relation between spheroid velocity dispersion and black hole mass. 

This is not the only problem for massive galaxies and their halos: there are others. For 
example, cooling flows occur in galaxy clusters at a rate that exceeds the observed rate 
of cool gas deposition. Additional entropy input, most likely in the form of preheating, 
is needed to reduce their role, in order to avoid an excessive rate of star formation in the 
brightest cluster galaxies. AGN are again the most likely culprit. Late feeding of AGN 
and jet-driven outflows provide a non-localised heat source that complements the early 
role of quasars in spheroid formation.. 

CDM-motivated theory does account for at least one observational result, namely the 
characteristic galactic mass of the Schechter function fit, L»~3x 10 10 L Q . This is based 
on the requirement that cooling within a dynamical time is a necessary condition for 
efficient star formation. The inferred mass-to-light ratio is reasonable, although one is 
left with the need to hide a significant fraction of the baryons. These are probably in the 
form of intracluster gas (Gonzalez, Zaritsky & ZabludofT 2007), within the observational 
errors, albeit some may be hidden at large radii. 

Modelling has succeeded in giving an approximate fit to the galaxy luminosity function. 
A relatively significant amount of dust is needed for the massive galaxies in addition to 
both quasar and AGN feedback, as well as reionisation and SNe feedback for the low 
mass galaxies (Bower et al. 2006; Croton et al. 2006; De Lucia & Blaizot 2007). The 
models are tuned to fit optical data. The fitting process is challenged as new, improved 
data becomes available. For example, the recent NIR data (Smith et al. 2008) reveals 
how fragile semi-analytical modelling has become. The UKIDSS SEDs yield improved 
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stellar mass estimates that are not particularly well fit by published models, for either 
low mass or massive galaxies. 

5. SEEKING A UNIFIED THEORY 

Feedback is essential for slowing down or even quenching star formation. One can 
make the case that feedback, in some circumstances, accelerates star formation. One of 
the most worrying issues is that of the IMF. Comparison of the UV rest frame-measured 
star formation history and the NIR-measured stellar mass assembly histories agrees at 
low redshift but diverges towards high redshift (Wilkins et al. 2008). The explanation is 
unlikely to be due to uncertainties in the dust extinction law at high redshift, where dust 
corrections are relatively small, A variation in the IMF is one possible explanation, the 
IMF becoming progressively top-heavy towards high redshift. 

There is a long history of top-heavy IMF discussions, largely centered around the G 
dwarf problem. This explanation is now largely discredited, in favour of a metal-free gas 
accretion history. The chemical abundance ratio signatures require that the Milky Way 
IMF was 'nearly invariant with time, place and metallicity (Gilmore & Wyse 2004). How- 
ever more recently, similar explanations invoking a top-heavy IMF have been advanced to 
explain the excess in the submillimctre galaxy counts (Lacey et al. 2008), the enrichment 
of the intracluster gas (Nagashima et al. 2005), the failure of semi-analytical modelling to 
reproduce the star formation efficiency enhancement towards high redshift (Dave 2008), 
the extragalactic light background (Fardal et al. 2007), and the luminosity and colour 
evolution of massive cluster galaxies (van Dokkum 2008). All of these interpretations, 
some of which have been given esoteric names by their proponents such as 'bottom-light' 
or 'paunchy', involve a partial suppression of the mass-carrying stars, below a solar mass, 
relative to those responsible for the bulk of the light and/or chemical yields. 

The best direct evidence for a truncated or top-heavy IMF was in the Arches star 
cluster (Stolte et al. 2005) but this disappeared with deeper data (Kim et al. 2007). Only 
within the central O.lpc of our galaxy, within the gravitational influence of the central 
supermassive black hole, is there a confirmed deficit of low mass stars relative to the 
observed massive stars (Figer 2008). This does perhaps support arguments for a top- 
heavy IMF in extreme situations, such as in the broad emission line regions of quasars 
where supersolar metallicities are observed and inferred to have been generated at z ^ 8 
(Mathias and Hamann 2008). However the gas masses involved, ^ 1O 4 M (Baldwin et 
al. 2003), seem to fall well short of what is needed if gas enriched by a top-heavy IMF 
were to have a global impact on galaxy evolution. 

Rather than add new parameters to the semi-analytical box of tricks, the future is 
more likely to bring greatly refined numerical treatments. It is clear that we are woefully 
resolution-limited when it comes to the 3-dimensional mixing of multiphase media that 
involves energy and momentum injection in one or more of the gas phases. Our under- 
standing of feedback is likely to change dramatically once our simulations match our 
aspirations. Consider the track-record: our insights into star formation, galaxy harass- 
ment and ram pressure stripping all underwent near 180 degree reversals once adequate 
adaptive grid power was brought to bear. I suspect that we are barely scratching the 
surface when it comes to understanding angular momentum transfer and AGN or even 
SNe feedback in a multiphase interstellar medium embedded within a protogalaxy. 
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